OPTIMAL CONTROL IN BOMBIERI'S AND 
TAMMI'S CONJECTURES 



Dmitri Prokhorov & Alexander Vasil'ev 

Abstract. Let S stand for the usual class of univalent regular functions in the unit disk U = 
{z : \z\ < 1} normalized by f(z) = z + a,2Z 2 + . . . in U, and let S M be its subclass defined by 
restricting < M in U, M > 1. We consider two classical problems: Bombieri's coefficient 

problem for the class S and the sharp estimate of the fourth coefficient of a function from S M . Using 
Lowner's parametric representation and the optimal control method we give exact initial Bombieri's 
numbers and derive a sharp constant Mo, such that for all M > Mo the Pick function gives the 
local maximum to |a4|. Numerical approximation is given. 



§1. Introduction 

Let S stand for the class of all holomorphic and univalent functions f(z) = z + Y^=2 a nZ n 
in the unit disk U = {z : \z\ < 1}. Its subclass of bounded maps \f(z)\ < M, M > 1, we denote 
by S M , S°° = S. During the long history of univalent functions the famous Koebe function 

oo 

K(z) = -^ = J2nz n eS 

^ n=l 

has been known to be extremal in many problems. A relevant sample is the most celebrated 
Bieberbach Conjecture [2] \a n \ < n proved by L. de Branges in 1984 [5,6]. In spite of many 
works about coefficient estimates in the class S, there are some difficult problems that are still 
unsolved, in particular, the Bombieri problem and the sharp upper bound for \a n \, n > 4, for 
the subclass S M that we will deal with. 

E. Bombieri [4] in 1967 posed the problem to find 

?i - Re a n 

a mn :=hminf , m,n>2, 

f^K m — He a m 

fes 
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/ — > K locally uniformly in U . We call a mn the Bombieri numbers. He conjectured that 
a mri = B mn , where 

n sin 9 — sm(n9) 
B mn = mm — — 

6>6[o,2ti-) m smv — sm{mt)) 

and proved that a mn < B mn for m = 3 and n odd. It is noteworthy that D. Bshouty and 
W. Hengartner [7] proved Bombieri's conjecture for functions from S having real coefficients in 
their Taylor's expansion. Continuing this contribution by D. Bshouty and W. Hengartner, the 
conjecture for the whole class S has been recently disproved by R. Greiner and O. Roth [9] for 
n = 2, to = 3, / e S. Actually, they have got the sharp Bombieri number a 32 = (e — l)/4e < 
1/4 = B 32 . 

It is easily seen that 043 = S43 = 023 = -B23 = 0. Applying Lowner's parametric representa- 
tion for univalent functions and the optimal control method we will find the exact Bombieri num- 
bers 042, cr 2 4, (T34 and their numerical approximations 042 ~ 0.050057..., 024 ~ 0.969556..., 
and (T34 ~ 0.791557... (the Bombieri conjecture for these permutations of to, n suggests 
B42 = 0.1, B 2 4 = 1, -B34 = 0.828427...). Of course, our method permits us to reprove the 
result of [9] about 032. 

Our next target is the fourth coefficient CI4 of a function from S M . An analogue of the Koebe 
function for this class is the Pick function 

00 

P M (z) = MK- 1 (K(z)/M) = z + J2Pn(M)z n . 

n=l 

The sharp estimate 1 0,2 j < 2(1 — 1/M) = p 2 (M) in the class S M is rather trivial and has 
been obtained by G. Pick [12] in 1917. The next coefficient 03 was estimated independently by 
A. C. Schaeffer and D. C. Spencer [14] in 1945 and O. Tammi [16] in 1953. The Pick function 
does not give the maximum to 1 0,3 1 and the estimate is much more difficult. M. Schiffer and 
O. Tammi [15] in 1965 found that |a 4 | < p 4 (M) for any f e S M with M > 300. This result was 
repeated by O. Tammi [18, page 210] in a weaker form (M > 700) and there it was conjectured 
that this constant could be decreased until 11. The case of function with real coefficients is 
simpler: the Pick function gives the maximum to 1 0.4 1 for M > 11 and this constant is sharp 
(see [17], [19, p. 163]). By our suggested method we will show that the Pick function locally 
maximizes 1 04 1 on S M if M > Mo = 22.9569 . . . and does not for 1 < M < Mq. This disproves 
Tammi's conjecture. 

§2. Preliminary statements 

The parametric representation of univalent functions is based on the Lowner differential 
equation and goes back to the famous Lowner's paper [11] where the author using an idea of 
semigroups of conformal maps derived the equation 



dw 
~dt 
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where the control function u = u{t) is piecewise continuous in t > 0. One finds the foundations 
of the parametric method, e.g., in [1,8,13]. It is convenient to make the change of variables 
t — > 1 — e~ l and rewrite Lowner's equation using the preceding notation for the independent 
variable as follows 

. dw —w e lu + w 

(1 » ¥ = (i-o^-. ' roU = *- 0£<£1 ' 

A. C. Schaeffer and D. C. Spencer [14] were the first who used Lowner's equation for the class 
S M and proved that the integrals 

(2) w = w(z,t) = (l-t)(z + a 2 (t)z 2 + ...) 
of the equation (1) represent a dense subclass of functions / G S M by 

(3) f(z) = Mw(z,l-l/M). 

Representation (3) is valid for all M > 1 including M = oo if the product in (3) is regarded as 
the limit as M — > oo. From now on, we will use the notation (2). 

We remark that the case u = 7r in (1) corresponds to the Koebe function in the class S or to 
the Pick function in the class S M by (3). Besides, the dense subclass of S M represented by (3) 
contains all functions that give the boundary points of the coefficient region 

V 4 M = {(a 2 , a 3 , Re a 4 ) : / G S M }, 1 < M < oo. 

For given real numbers fi and u, we will consider the linear functional 

L(/i, u; /) = a 2 + /ua 3 + ua 4 

in S M . The Koebe function K(z) maximizes Re L(0, 0; /) in S, and similarly, the Pick function 
Pm(z) maximizes Re L(0, 0; /) in S M . We will describe the set of \i and v for which the local 
maximum of Re v\ f) in S M is attained by the Pick function and will apply this result to 
the extremal problems for the class S or S M stated in the introduction. 

We write Ofe(t) = X2ks(t) +ix2k-2(t), k = 2,3,4. Substituting (2) into (1) we obtain the 
following differential equations 

x\(t) = — 2cos-u, xi(0)=0, 
&2(t) = 2 sin u, 0:2(0) = 0, 

(4) xs(t) = — 4(xi cos u + x 2 sinw) + 2(t — 1) cos2w, x 3 (0) = 0, 
x±{t) = 4(xi sin-u — x 2 cos u) — 2(t — 1) sin2-u, ^4(0) = 0, 
xs(t) = — 2((2x 3 + x\ — x 2 ) cos u + 2(x4 + X\x 2 ) sin u) 

+ 6(t - l)(xi cos 2u + x 2 sin 2u) - 2(t - l) 2 cos3-u, x 5 (0) = 0. 
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The extremal problem 

Re L(fi, v\ f) — > max 
in the class S M is equivalent to the extremal problem 

(5) xi (1 - 1/M) + nx 3 (l - 1/M) + vx 5 (l - 1/M) -> max 

for solutions to the system (4). The parametric representation (4) for the coefficients generated 
by Lowner's equation allows us to apply the classical variational methods [3] or Pontryagin's 
maximum principle [10]. We introduce the Hamiltonian function in order to formulate the 
necessary extremum conditions for the problem (5) 

H(t,x, $,u) = -2 cos u^i + 2 sin-u^ 

— (4(xi cos u + x 2 sin it) — 2(t — 1) cos2w)\I/3 

(6) + (4(xi sin u — X2 cos it) — 2(t — 1) sin 2u)^4 

— (2 ((2x3 + x\ — x\) cosu + 2(x4 + X1X2) sin it) 

— 6(t - l)(xi cos 2it + x 2 sin 2it) + 2(* - l) 2 cos 3u)W 5 , 

where x = (xi, . . . , a?5) T satisfies (4) and \& = (\I/i, . . . , ^$) T satisfies the conjugate system 

^1 = Acosu^s — 4 smu^ 4 + (4xi cos it + 4x 2 sin it — 6(t — 1) cos2it)* 5 , 
4> 2 = 4sinw\E'3 + 4 cos it \I/ 4 — (4x 2 cos-u — 4xi sin it + 6(t — 1) sm2u)^ 5 , 

(7) #3 = 4cosit# 5 , 
^4 = 4sin-u^ 5 , 
^5 = 0, 

and the transversality conditions 

*i(l - 1/M) = 1, * 3 (1 - 1/M) = (jl, #5(1 - 1/M) = v, 
^ > * 2 (1 - 1/M) = #4(1 - 1/M) = 0. 

The optimal control function u* corresponding to the extremal function /* in (5) satisfies 
Pontryagin's maximum principle 

(9) max#(t,£*,^*,w) = H(t,x*,V*,u*), < t < 1 - 1/M, 

u 

where (x*, ^*) is the solution to (4) and (7) with it = it* in their right-hand sides. Hence, it* is 
a root of the equation 

(10) H u (t,x,V,u) = 
for x = x* and * = 
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Lemma 1. Let us suppose that a control function u in (4) and (7-8) generates the solutions 
x(t) and ^(t) for which u satisfies (9), is unique up to the 2ir -translation, and 

(11) H uu (t,x,V,Tr)^0, 0<t<l-l/M. 

Let us denote by (x(t, £), ^(t, £)) solutions to (4) and (7) with the initial conditions \I/(0, £) = 
^(0) + £ and u = it(t, £) in their right-hand sides satisfying the maximum principle (9). Then 
for £ — > 0, we have the following asymptotic behaviour 

|| (x(l - 1/M, 0, *(1 - 1/M, 0) - (x(l - 1/M), *(1 - l/M))\\ = o(l), 

where \\ ■ || is the Euclidean vector norm. 

Proof. Since there exists a unique solution u satisfying (9) and (11), the same is true for a slightly 
changed parameters of the function H. Therefore, equations (9) and (10) locally determine a 
unique continuous implicit function u = u(t, x, \1/) satisfying the maximum principle. Writing 
u(t,£) = u(t, x(t, £), ty(t, £)) we substitute it into (4) and (7). Now we apply the theorem on 
the continuous dependence of solutions of differential equations on the initial conditions and 
complete the proof of Lemma 1. □ 

If the Pick function Pm is extremal for (5), then u = n is the optimal control function, (4) 
and (7) give (x(t),V(t)) = (x°(t),V°(t)), where 

(12) x° 1 (t)=2t, x° 3 (t) = 5t 2 - 2t, x%(t) = x1(t) = 0, 
and 

"( t - 1 + s) 2 + (^- 8 "- 4 ")( t - 1 + ^) + 1 ' 

-4v(t-l + ^+p, *§(<)= v, *»(t) = $»(() = 0. 

The conditions of Lemma 1 play the key role as a necessary local extremum condition for the 
Pick function Pm- To verify these we substitute (12) and (13) into the Hamiltonian function 
H(t, x, u) given by (6) and study the extremum properties of H(t, u) = H(t, x°(t), ^°(t), u) 
which is just a cubic polynomial of cos u. Let us describe a set of suitable real parameters (//, v) 
satisfying Lemma 1. Let D(M) denote the maximal domain in the (fx, z/)-plane which is starlike 
with the respect to the origin and satisfies the following conditions: 

[i] H(t, u) as a function of y = cos u attains its maximum on [—1, 1] only at y = — 1 for all 
te [0,1 -1/M]; 

[ii] H uu (t,7r)^0, 0<t<l-l/M. 

We will consider (/i, v) e D(M). The point x°(l - 1/M) belongs to the boundary <9V 4 M of V 4 M 
and is given by the Pick function. Each x E dVl 4 can be obtained as a solution to the system 
(4) with a certain optimal control function u. Lemma 1 admits a reverse formulation. 



(13) 



*?(*) 
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Lemma 2. Let (n,u) e D(M) and denote by (x~(t, £), \I> , £)) solutions to (4) and (7) with 
the boundary conditions — 1/M, £) = \1/ (1 — 1/M) + £. Let it = u~(t, £) m t/iezr right-hand 
sides satisfy the maximum principle (9). If 

||z-(l-l/M,O-z (l-l/M)|| =o(l), as^O, 

thcfi 

||vl>-(0,0-f (0)||=o(l). 

The proof of Lemma 2 is similar to that of Lemma 1 reversing the direction of variation of t 
from 1 - 1/M to and noting that x~(0) = x°(0) = 0. 

Lemmas 1 and 2 imply that if (x(t),^(t)) is given by (4) and (7), (x(l - 1/M), #(1 - 1/M)) 
is close to (x°(l - 1/M), *°(1 - 1/M)), and x(l - 1/M) e dV 4 M , then (x(t), #(£)) is equal to 
(x(t, £), £)) for a certain £ close to 0. 

The principles of calculus of variations interpret geometrically the transversality conditions 
as an orthogonality property of \1/(1 — 1/M) to all possible variations of x(l — 1/M) in dV^ 1 . 
Nevertheless, we rigorously prove this fact for completeness. 

Lemma 3. Let us suppose that (fi, v) G D(M) and the initial conditions in (7) are 

(14) v|/(0,O = * O (0)+£, £ = ee, e > 0, e = (e u . . . , e 5 ) T , ||e|| = 1, 
and 

(15) x{t,£) = x°(t) + e8x(t) + o(e), e -> 0. 
T/ien *°(1 - 1/M) is orthogonal to Sx(l - 1/M), if 8x(l - 1/M) ^ 0. 

Proof. First we note that the conditions of Lemma 1 guarantee the differentiability of x(t,ee) 
with respect to e at £ = 0. Therefore, the representation (15) is valid and according to Lemmas 1 
and 2 the expansion (15) produces all possible variations Sx(l — 1/M) at x°(l — 1/M) associated 
with #°(1 - 1/M). 

We denote the column of the functions in the right-hand side of (4) by g(t,x, u) and rewrite 
(4) in the vector form as 

(16) x = g(t,x,u), x(0)=0. 
The system (7) is equivalent to 

(17) i? = - — (t,x,V,u), tf(0) = V°(0)+ee T . 

ox 

Substituting (15) into (16) we obtain 

dSx 

— — = g x dx + g u u £ , 



OPTIMAL CONTROL. 



7 



by differentiating with respect to e at e = 0. This, together with (4), (6), (10), and (17), imply 
that 

j t ({5x) T V) = (6x) T gT* + u e gTy + (dxfv = 0, 

because the second term is equal to u £ H u (t, x, ^, u) = 0, and the remaining terms give the zero 
sum since 4/ — —g^^. 

Thus, we see that (5x) T ^ does not depend on t and vanishes at t = because 5x(0) = 0. 
This completes the proof of Lemma 3. □ 

§3. Local extremum conditions 

Lemma 4. Under the conditions of Lemma 3 we suppose that in the vector e = (ei, . . . , e^) T , 
which corresponds to the variation of ^°(0) in (14), the coordinates e 2 and e^ vanish. Then 
Sx(l - l/M) = 0. 

Proof. The condition e 2 = = implies that the systems (4) and (7) have vanishing coordinate 
solutions x 2 (t) = x 4 (t) = and ^ 2 {t) = ^4^) = 0. In this case the Hamilton function if is a 
polynomial of y = cos-u, which has a unique maximum on [—1, 1] at y = —1. Its derivative with 
respect to y does not vanish at y = — 1 for e > sufficiently small. Hence, u = n is a unique 
optimal control function for such e and x(l — l/M, £) = x°(l — l/M), that ends the proof of 
Lemma 4. □ 

By analogy with the expansion (15) in Lemma 3 we have the expansion 

Lemma 4 shows that the condition e 2 = = implies 5x = and 5^> 2 = <^4 = 0. Only 
\I/3, and ^5 can vary in this case . It follows from Lemmas 1-4 that we should consider 
variations by (14) with e\ = e^, = e$ = in order to study the character of the point 

- l/M) edV 4 



M 



4 



Let us set £ = (0, p, 0, q, 0) T with arbitrary real p and q, and study x(l — l/M, £) in a 
neighborhood of x°(l — l/M). In other words, we will solve the systems (4) and (7) with the 
initial conditions x(0) = and (14), which we rewrite with coordinates 



M 



Let 



*i(0) 


= 3i/ ^1 


* 2 (0) 


= P, 


* 3 (0) 


= 4u (l 


*4(0) 


= Q, 


* 5 (0) 


= v. 


: (p, g) 





xi(l - l/M) + //x 3 (l - l/M) + i/x 5 (l - l/M) 
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be a real valued mapping from the (p, (j)-plane onto the linear combination of the components of 
the solution to the Cauchy problem for the systems (4) and (7) with the initial conditions (18). 
The control function u in the right-hand side of (4) and (7) satisfies the maximum principle. 
The mapping F is well defined in a neighborhood of (0,0) if (//, v) G D(M). In this case 
u = u(t,x^) is an implicit function defined by (10). Since (x,^) depends only on (p,q), we 
denote by it(£, p, q) = u(t, x(p, q), ^{p, q)). We note that F(0, 0) = Re L(/x, u; Pm) and the values 
of F{p, q) correspond to those of Re L{p>, v\ f) with respect to the variations of Pm generated 
by the initial conditions ^(O) = p and ^4(0) = q. 

Theorem 5. Let us suppose that (/J,,v) G D(M). If Pm locally maximizes Re L(fi,i/; /) in 
S(M), then 

F pp (0, 0) < 0, F pp (0, 0)F qq (0, 0) - F p 2 g (0, 0) > 0. 

Conversely, if 

F pp (0, 0) < 0, F pp (0, 0)F qq (0, 0) - F p 2 g (0, 0) > 0, 
then Pm locally maximizes Re L(/j,, v\ f) in S(M). 

Proof. We first claim that (0, 0) is a critical point of F(p, q). Indeed, substituting u = u(t,p, q) 
in the three equations in (4) (for x\, X3, x§) and differentiating them with respect to p and q, we 
obtain differential equations for (xk) P and (xk) q , k = 1,3,5, with vanishing initial conditions. 
Substituting there p = q = and u = n, X2(t) = x±(t) = we find that all derivatives (xk) p and 
(xk) q , k = 1, 3, 5 are identically zeros and, hence, F p (0, 0) = ^(0, 0) = 0. 

The first statement of Theorem 5 means that the quadratic form of the second differential of F 
at (0, 0) is negatively semi-definite which is the necessary condition of local extremum. Similarly, 
the second statement signifies that the above quadratic form is negative definite which is the 
sufficient local extremum condition. This completes the proof of Theorem 5. □ 

The same reasoning can be made for the class Sf^ of functions / G S M with real Taylor 
coefficients a n , n > 2. The coefficients of an arbitrary boundary function fn G for the set 
V^l = {(02,03,04) : / G Sj^} can be obtained by integrating the systems (4) and (7) with 
the control function u satisfying the maximum principle (9) and with vanishing initial values of 
\l/2(0) and ^4(0). Therefore, variations of ^2(0) and ^4(0) are forbidden and the Pick function 
P M e locally maximizes L(fi, u; f) in if ( /U , u ) G D{M). 

Now we will apply Theorem 5 to construct an analytic and numerical solution process. We 
need to calculate the partial derivatives u p and u q at (0, 0) to evaluate the partial derivatives of 
F at (0,0). Differentiating (10) with respect to p and q we obtain 

HuxXp + Hu^f^p + H uu u p = 0, 

HuxXq + H u ^ q + H uu U q = 0, 

which leads us to the formulae 
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and 

H ux x q + H u ^ t 



(20) u q =- — 



ltu f - 4 I 2, + t ~ »)* + *> + 1 - 4(2 "J + 



Direct calculation gives 

(21) H m (t,A*°,r) = -2 . M ,j M - M2 

Differentiating (6) with respect to corresponding variables at u = ir we find 

(22) H UX2 (t, x°, vl/ , tt) = 4 (V (/ + f - + n 

(23) H UX4 (t,x°,-$°,ir) =4i/, 

(24) ^ 2 (t,x°,vp°,7r) = -2, 

(25) ^ 4 (t,x°,vp°,7r) =4(1-3*). 

Differentiating (7) with respect to p and g at (0, 0) with u = ir, (xi) p (t) = {x\) q {t) = x 2 (t) = 
we see that 

(26) (*!)„(*) = = (*s) P (f) = (*s),(f) = (*5)p(«) = (*5) fl («) = 0. 

The formulae (22-26) allow us to calculate the numerators in (19) and (20) as 

4 ^ (t + 1 - + f?j (x 2 ) p + 4z/(x 4 ) p - 2(tf 2 ) p + 4(1 - 3*)(tf 4 ) P 

and 

/ / 4 \ \ , N 



respectively. 

From (4) and (7) we conclude that \l/ 4 = 2ux2 yields the equalities (^a) p = 2z/(i 2 ) p and 

(\j> 4 ) = 2v(x 2 ) q . The initial conditions (* 4 ) p (0) = and (* 4 ) g (0) = 1 imply that (* 4 ) p = 

2z/(x 2 ) p and (* 4 )q = 2v(x2) q + 1. 

We substitute the last relations and (21) into (19-20), and finally, get 

( (3 - 5£ - 4/M) 1/ + fj.) 2y 4 + 2^5 - ye 

(27) -Up = 



16ut 2 - 4(2z/ + 4i//M - fj)t + 2v + 1 - 4(2i/ + //)/M + 15z//M 2 
where y 4 := (x 2 ) p , 2/5 := Mp, 2/6 := (*2)p, and 

((3 - 5/ - A/M)v + fi)2y 10 + 2uy n - y 12 + 2(1 - 3/) 



(28) u q = 



16z4 2 - 4{2v + Au/M - n)t + 2v + 1 - 4(2i/ + //)/M + 15^/M 2 ' 
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where y 10 := (x 2 ) q , 2/n := (x 4 ) g , 2/12 := (*2)g- 

Set yi := (xi) pp , y 2 := (x 3 ) pp , and 1/3 := {x$) pp . Differentiating (4) twice with respect to p 
at u = 7r, :r 2 = 2:4 = \I/ 2 = ^4 = (^i)p = (^3)? = we obtain 

(29) yi = -2u 2 p , yi(0) = 0, 

(30) y 2 = 4(y 1 + 2y 4 « P -2(2t-l) U J), 2/ 2 (0) = 0, 

2/ 3 = 2(7* - 3)yi + 4y 2 - 4y| + 8(5* - 3)y 4 w p + 8y 5 u P 

(31) -2(47t 2 -46t + 9)wJ, 2/3 (0) = 0. 

Similarly, we differentiate the remaining equations in (4) with respect to p and obtain 

(32) y 4 = -2u p , 2/ 4 (0) = 0, 

(33) y 5 = 4(y 4 + (l-3t)u p ), y 5 (0) = 0. 

Finally, we differentiate the second equation in (7) with respect to p and get 



■h = -4" (j(t + 1 ~ + I*/") % + 2/4^ , 2/e(0) = 1. 



(34) y 6 

Summarizing, we have deduced an evaluation algorithm for F pp (0, 0) expressed by the following 
theorem. 

Theorem 6. Suppose (fx, v) G D(M). Let yi(t), . . . , ye{t), < t < 1 — 1/M, be solutions to the 
Cauchy problem for the differential equations (29-34)- Then the relation 

F pp (0, 0) = - 1/M) + fjiy 2 (l - 1/M) + uy 3 (l - 1/M) 

is valid. 

Remark. The subsystem (32-34) can be solved independently because these equations do not 
contain 2/1, y 2 , and y 3 . 

Calculation of F qq (0, 0) and F pq (0, 0) may be handled in much the same way. Let 2/7 := (xi) qq , 
y 8 := (x 3 ) qq , and y 9 := (x 5 ) qq . From (4) we have 

(35) y 7 = -2u 2 q , y 7 (0) = 0, 

(36) y 8 = 4(y 7 + 2y 10 u q -2(2t-l)u 2 q ), y 8 (0) = 0, 
y 9 = 2(7* - 3)y 7 + 4y 8 - 4y 2 10 

(37) + 8(5* - 3)y 10 u q + 8 Vll u q - 2(47* 2 - 46* + 9)u 2 , y 9 (0) = 0. 

Differentiating the two even equations in (4) and the second equation in (7) with respect to q 
we get 

(38) y w = -2u q , yio(0) = 0, 

(39) S/ii = 4(y 10 + (l-3*K), 2/11 (0) = 0, 

(40) y 12 = -4v(U+l-^ + fi)u q + yi o)-4, 2/12(0) = 0. 
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Let 1/13 := (xi) pq , J/14 := {x^) pq , and j/15 := {x^) pq . We continue in this fashion differentiating 
(4) with respect to p, and subsequently, with respect to q and obtain 

(41) m = -2u p u q , y 13 (0) = 0, 

(42) yi4 = 4(yi 3 + y 4 u q + yi u p ) - 8(2t - l)u p u q , 2/14(0) = 0, 
y 15 = 2(7t - 3)2/13 + 4y 14 - 4z/ 4 z/io + 4(5* - 3)(z/ 4 w g + yi « P ) 

(43) + 4y 5 u q +4y llUp -2(47t 2 -46t + 9)u p u q , y 15 {0) = 0. 

Summing up the calculation process we derive the following theorem. 

Theorem 7. Suppose (/U, z/) G D(M). Let yr(t), . . . ,yi2(t), < t < 1 — 1/M, be solutions to 
the Cauchy problem for the differential equations (35-40)- Then, the relation 

F qq (0, 0) = y 7 (l - 1/M) + m (l - 1/M) + uy 9 (l - 1/M) 

holds. Let y/t(t), y$(t), ye(t), and yi3(t), yi4(t), yis(t), < t < 1 — 1/M, be solutions to the 
Cauchy problem for the differential equations (32-34) an d (41~4$) respectively. Then, the rela- 
tion 

F pq (0, 0) = y 13 (l - 1/M) + ^4(1 - 1/M) + uy 15 (l - 1/M) 

holds. 

Remark. As in the remark after Theorem 6 we note that the subsystem (38-40) can be solved 
independently because these equations do not contain 2/7, y 8 , and yg. 



§4. Explicit integration 

An explicit integration of the systems in Theorems 6 and 7 is possible only in the case v = 0, 
i.e., when the linear functional L does not depend on CI4. In this case the two last equations in 
(4) and (7) disappear and the mapping F becomes a function of a variable p as 



F : p -> xi(l - 1/M) + fjLx 3 (l - 1/M). 

The criterion in Theorem 5 is reduced to the inequality F"(0) < 0, where F"(0) = yi(l — 1/M) + 
^2/2(1 — 1/M). The systems (29-34), (35-40), and (41-43) are reduced to the four equations 
(29), (30), (32), and (34). We substitute v = into (34, 27) and obtain that 



(44) 
and 
(45) 



y 6 = -4iiu p , y 6 (0) = 1, 



U n = 



2^2/4 - 2/6 

4fit + 1 -4/x/M' 
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Compare (32, 44) and observe that y§ = 1\xy^ which implies 

This equation allows us to exclude y@ from (45) and integrate the system of equations (29), (30), 
and (32) with 

(46) u p - - — 



Aiit + 1 - Ali/M 
The equations (29) and (32) give 

^ Vl ® = Yli (l -4///M + 4//* ~ 1-All/m) 

and 

1 , 1 - Au/M + Aut 

48 yAt) = — log ^ ; — . 

V 1 y W 2// & 1-4///M 

Substituting (46) and (48) into (30) we integrate it, then, taking into account (47), we finally 
obtain that 



(49) yi (t) + fj,y 2 (t) = 



log 2 ((l + 4/x - 8/i/M)(l - Atx/M)) + log 1 + !// "^ //!/ 



1 - 4///M 

Making use of (49) we formulate the following theorem as a corollary of Theorem 6. 
Theorem 8. Suppose (/i, 0) G D(M) and li satisfies the inequality 



(50) — 
V ; 2// 



log 2 ((l + 4/x - 8/*/M)(l - 4/VM)) + log 1 V- 4 



> 0. 



Then, the Pick function Pm locally maximizes Re L(/x, 0; /) = Re (a2 + ^03) in S M . If the 
left-hand side of (50) is negative, then Pm does not give a local maximum to Re L(/i,0;/) in 
S M . 

We note that {li : (/z, 0) G D = D(oo)} = {li : \i > —0.25}. If M = 00, then the inequality 
(50) is of the form 

^-(1 + log(l + All)) log(l + 4/z) < 0, 

that is equivalent to the best possible inequality /U > — Ao := — (e — l)/4e. This means that the 
Bombieri number CT32 is equal to Ao- This result has been obtained recently by R. Greiner and 
O. Roth [9]. 
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§5. BOMBIERl'S NUMBER <7 42 

Now we apply Theorems 5-7 to evaluate Bombieri's number 042. 
Proposition. Let m, n > 2 be fixed integers. Then 

o~mn = sup{A G K : Re (a n — Aa m ), is locally maximized on S by K(z)}. 

The proof of this statement is quite obvious and can be found, e.g., in [9]. 
According to Theorem 5, Bombieri's number 042 is calculated as 

- inf{z/ : F pp (0, 0) < and F pp (0, 0)F qq (0, 0) - F p 2 g (0, 0) > for // = 0, M = 00, v e [1/, 0]}. 

Theorems 6 and 7 reduce the problem to the solution of the equations 

(51) yi(l) + ^ 3 (l) = 0, 

or 

(52) + i/y 3 (l))(y 7 (l) + ^9(1)) - (2/13(1) + ^i 5 (l)) 2 = 0, 

where yi (£),... , yis(t) are solutions to the Cauchy problem for the differential equations (29-43) 
with /j = and M = 00. Thus we are able to formulate the following result. 

Theorem 9. Bombieri's number 042 is equal to the maximum of the negative roots of the 
equations (51) and (52) multiplied by (-1), where yi(t), . . . , yis(t) are solutions to (29-43) for 
[i = and M = 00. 

To illustrate this result we give the numerical approximation 042 ~ 0.050057... 

Remark. We used Wolfram's Mathematica to evaluate numerically 042 as well as other Bombieri's 
numbers in the next sections. More precisely, the combination of 4-th order Runge-Kutta and 
Adams methods is used. It allows to reach a higher precision at short computational time. 
Actually, the level of precision imply the number of iterations and can be prescribed as the 
machine precision. For solving the above systems the level 10~ 12 is reached in time of order of 
few minutes on a machine with a Pentium 4 processor (1800 Mhz) and having 512 Mb of RAM. 

Remark. Note that \y : (0,v) e D = D(oo)} = {v : v > —0.1}. O. Roth showed (private 
communication) that 042 < 0.050284..., the latter number is equal to (2 — Re a2)/(4 — Re 04) 
for a function sequence which is critical in the problem of finding 0*32- 
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§6. BOMBIERl'S NUMBERS (J 2 4 AND (J34 

The problem of finding Bombieri's number (T24 is reduced to description of v for which the 
Koebe function K minimizes Re L(0, u; f) in S. To follow the preceding scheme we will write 
L 1 (z/; /) = 04 + z/a2 in place of L. Now cr 24 is equal to the supremum over all real values A24, 
such that Re L x (— A24; /) is locally maximized by the Koebe function K in S. Again we consider 
the system of equations (4), the Hamilton function H(t, x, u) given by (6), and the system of 
equations (7), where the transversality conditions (8) are replaced by 

(53) *i(l) = ^, *s(l) = l, # 2 (1) = * 3 (1) = *4(1) = 0. 

For u = 7T let us denote the integrals for (4) and (7) by (x°(t), where x\(t), . . . , x^(t) are 

given by (12) and 

(54) q{(t)=t 2 -10t + 9 + v, = -4(t-l), = 1, *2(*) = *4(*) = 0- 

As before, the control function u = u(t) satisfies the maximum principle (9) and, hence, the 
equation (10). We consider only v e D 1 , where D 1 denotes a maximal interval on the z/-axis 
that satisfies the following conditions: 

[i] H x {t, u) = H(t, x°(t), \& 1 (t) , u) as a function of y = cos u attains its maximum on [—1, 1] 
only at y = — 1 for all t G [0, 1]; 

[ii] Hl u (t,7r)^0, < t < 1. 

Note that D 1 = {z/ : z/ > -1}. 

Similarly to Section 3 we solve the systems (4) and (7) with the initial conditions x(0) = 
and 

(55) *!(0) = 9 + ^ * 2 (0)=p, * 3 (0) = 4, * 4 (0) = g, * 5 (0) = 1, 

where p and q are arbitrary real numbers. The initial conditions (55) are thought of as variations 
of (54) at t = 0. 
Let 

F 1 : {p,q)^v Xl {l) + x 3 {l) 

be a real valued mapping from the (p, </)-plane to the linear combination of the components 
of a solution to the Cauchy problem for (4) and (7) with the initial conditions (55) and the 
control u satisfying the maximum principle. The mapping F 1 is well defined in a neighborhood 
of (0,0) if v e D , and u = u(t,x,^>) is an implicit function defined by (10). We denote by 
U (^P7<?) = u (ti x {Pi q)i ^{Pi <?))• Lemmas 1-4 and Theorem 5 admit an analogous formulation 
for the functional L x {v; f). So, if K locally maximizes Re L 1 (z/; /) in S for v e D 1 , then 

F^(0, 0) < 0, F^(0, 0)F^(0, 0) - (Fi) 2 (0, 0) > 0. 

Conversely, if 

Fi p (0, 0) < 0, Fi p (0, 0)F* (0, 0) - (F„M 2 (0, 0) > 0, 
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then K locally maximizes Re L 1 (z/; /) in S. Direct calculation gives 

(56) H uu (t, x°, #\ tt) = -2(16t 2 - 8* + v + 2), 

(57) ^ X2 (t,x°,^ 1 ,7r)=4(t + 1), H UX4 (t,x°,V 1 ,n)=4, 

(58) ^ 2 (t,x°,* 1 , 7 r) = -2, ^ 4 (t,a; ,* 1 ,7r) = 4(l-3t). 

The formula (26) remains true as well as (^4)^ = 2(x 2 ) p and (^4)3 = 2{x 2 ) q + 1. Substituting 
the last relations and (56-58) into (19) and (20) we obtain 

/con ( 3 - 5t)2y 4 + 2y 5 - y 6 

(59) -Up = 



(60) u 



m 2 - 8t + v + 2 ' 
_ (3 - 5t)2y 10 + 2yn - y 12 + 2(1 - 3t) 
g 16t 2 - 8t + v + 2 



Evidently, the formulae (29-33), (35-39), and (41-43) are valid for our case, whereas the equa- 
tions (34) and (40) are transformed into 

(61) y 6 = -4(t + l)up - 4y 4 , ye (0) = 1, 

(62) y 12 = -4(t+l)u q -4y w -4, y 12 (0) = 0, 

respectively. 

Summing up above calculation for the evaluation algorithm we state that by analogy with 
Theorems 6 and 7 the problem of finding 024 is reduced to the solution of the equations 

(63) ^ 1 (l) + y 3 (l) = 
or 

(64) (u yi (l) + y 3 (l))(isy 7 (l) + y 9 (l)) - (uy 13 (l) + y 15 (l)) 2 = 

where yi(t), . . .,yis(t) are solutions to the Cauchy problem for the differential equations (29- 
33), (61), (35-39), (62), and (41-43) with u p and u q given by (59) and (60). We formulate the 
following result. 

Theorem 10. Bombieri's number 024 is equal to the maximum of negative roots of the equations 
(63) and (64) multiplied by (-1), where yi(t), ... ,y 15 (t) are solutions to (29-33), (61), (35-39), 
(62), (41-43), and (59-60). 

Similarly to Bombieri's number 042, the numerical approximation for o~ 2 4 is 0.969556..., which 
is the maximal negative root of (64) . 

In order to evaluate 034 we consider the functional N(fi; /) = 04 + /UCI3 as soon as 034 is 
equal to the supremum over all real values A34 for which Re iV(— A34; /) is locally maximized by 
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the Koebe function K in S. Now the systems (4) and (7) are supplied with the transversality 
conditions 

(65) * 3 (1) = ^, * 5 (1) = 1, *i(l) = * 2 (l) = *4(l) = 0. 

For u = it we denote the integrals of (4) and (7) by (x°(t), \I/ 2 (t)), where Xi(t), . . . , x\{t) are 
given by (12) and 

V 2 (t)=t 2 - (10 + 4/^ + 9 + 4^, *i(t) = -4(*-l)+//, 

*§(*) = i, = = o. 

We consider only G -D 2 , where D 2 is a maximal interval on the jU-axis which satisfies the 
conditions: 

[i] H 2 (t : u) = H(t, x°(t), \E f2 (t), u), as a function of y = cos-u, attains its maximum in [—1, 1] 
only at y = — 1 for all £ G [0, 1]; 

[ii] U 2 Jt,7T)^0, < t < 1. 

Note that D 2 = {/j : fi > -2( v / 2 - 1)}. 

Again we solve the systems (4) and (7) with x(0) = and 

(67) *i(0) = 9 + 4//, * 2 (0)=p, # 3 (0) = 4 + ^, * 4 (0) = g, * 5 (0) = 1, 

where p and q are arbitrary real numbers. The conditions (67) are thought of as variations of 

(66) at t = 0. 
Let 

F 2 : (p,q) -> nx 2 (l) + x 3 (l) 

be a real valued mapping from (p, q) -plane into the linear combination of components of the 
solution to the Cauchy problem for (4) and (7) with the initial conditions (67) and u satisfying the 
maximum principle. We preserve preceding denotations for u(t,p, q), and similarly to Theorem 
5, for fx G -D 2 , we assert, that if K locally maximizes Re N(fi; f) in the class S, then 

F%% 0) < 0, F 2 p (0, 0)F 2 ? (0, 0) - (F 2 g ) 2 (0, 0) > 0. 

Conversely, if 

F 2 p (0, 0) < 0, F 2 p (0, 0)F 2 g (0, 0) - (F 2 g ) 2 (0, 0) > 0, 

then K locally maximizes Re N(fj,; f) in S. 
Direct calculation gives 

(68) H uu {t, x°, tf 2 , tt) = -4(8t 2 - (4 - 2fi)t + 1), 

(69) H UX2 (t,x°^ 2 ,ir) = 4(t + l + iu), H UX4 {t, x°, tf 2 , tt) = 4, 

(70) H u v a (t,x ,* 2 ,ir) = -2, i^ 4 (t,x°,tf 2 ,7r) =4(1-3*). 
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The formulae (19), (20), and (68-70) lead to 

(3 - 5t + /j)2y 4 + 2y 5 - y 6 



(71) u p = 

(72) u q 



IQt 2 - (8 - 4/i)t + 2 
(3 - 5t + /i)2y 10 + 2y n - y 12 + 2(1 - St) 



16t 2 - (8 - 4/j)t + 2 
The formulae (34) and (40) are changed to 

(73) y 6 = -4(£ + 1 + /jl)u p - 4y 4 , y 6 (0) = 1, 

(74) y 12 = -4(t + 1 + fi)u q - 4y 10 - 4, y 12 (0) = 0, 

respectively. 

Let us sum up the results and state, that similarly to Theorems 6 and 7, the problem of 
finding 034 is reduced to the solution of the equations 

(75) ^ 2 (l) + y 3 (l) = 0, 
or 

(76) (/^(l) + y 3 (l))(^ 8 (l) + 2/9(1)) - ( m4 (l) + t/ 15 (l)) 2 = 0, 

where yi(t), . . . , 2/15 (t) are solutions to the Cauchy problem for the differential equations (29-33), 

(73) , (35-39), (74), and (41-43) with u p and u q given by (71) and (72). 

Theorem 11. Bombieri's number 034 is equal to the maximum of the negative roots of equations 
(75) and (76) multiplied by (-1), where y^t),..., y 15 (t) are solutions of (29-33), (73), (35-39), 

(74) , (41-43), and (71-72). 

Numerical methods applied to (75) and (76) show that 034 w 0.791557... 

§7. Fourth coefficient of bounded univalent functions 

In this section we will find a sharp constant Mo, such that for all M > Mo the Pick function 
Pm gives the local maximum to Re a 4 (and to | a-4 1 ) over all univalent functions / G S M . In the 
class Sft C S M of functions with real Taylor coefficients the Pick function gives the maximum 
to CI4 for M > 11 and this constant is sharp (see [17], [19, p. 163]). That is why we can consider 
only M > 11. For this problem an analog of Theorem 5 holds and allows us to construct an 
evaluation algorithm. Therefore, we again consider the systems (4) and (7) supplied now with 
the transversality conditions 



(77) * 5 (1-1/M) = 1, #i(1 - l/M) = ■■■ = #4(1 - l/M) = 0. 
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For u = 7r, let us denote the integrals of (4) and (7) by (x°(t) : \E f3 (t)), where 

*f (t) = t 2_ A Q _ 16V g _ 24 15 ^3 (t) = _ 4 L_ l 1 

(78) V M / M M 2 V M 

*i(t) = i, *i(t) = *ut) = o. 

Varying conditions (78) at t = we put 

(79) * 1 (0) = 9-^ + ^, * 2 (0)=p, *3(0) = 4^1-^V ¥ 4 (0) = ?J * 5 (0) = 1, 

where (p, </) G M 2 . 
Let 

F 3 :(p,Q)-^a;5(l-l/M) 

be a mapping from the (p, g)-plane into the axis of the third component of the solution to the 
Cauchy problem for (4) and (7) with the initial conditions (79) and u satisfying the maximum 
principle. Preserving notations of Section 6 for u(t,p,q), similarly to Theorem 5, we observe 
that for M > 11, if Pm locally maximizes Re 04 in S(M), then 

F p 3 p (0, 0) < 0, F p 3 p (0, 0)F 3 g (0, 0) - (F p 3 g ) 2 (0, 0) > 0, 

and conversely, if 

F p 3 p (0, 0) < 0, F p 3 p (0, 0)F 3 g (0, 0) - (F p 3 g ) 2 (0, 0) > 0, 

then Pm locally maximizes Re 04 in S(M). 
Direct calculation gives 

(80) H uu (t, x°, vP 3 , n) = -2 (V - (s + t + 2 - ± + J* ) , 

(81) H UX2 (t,x°^ 3 ,ir)=4(t + l-j^y H ux ,(t,x°,^.7r) = 4, 

(82) if^ 2 (t,x°,* 3 ,7r) = -2, H u < Su (t,x ,V 3 ,Tr) = 4(l-3t). 

The formulae (19), (20), and (80-82) lead to 

(3-4/M -5t)2y 4 + 2y 5 - y 6 



(83) u p = 

(84) u q = 



16t 2 - (8 + 16/M)t + 2- 8/M + 15/M 2 ' 
(3 - 4/M - 5t)2y 10 + 2yu - m + 2(1 - 3g 
16t 2 - (8 + 16/ M)t + 2 -8/M + 15/M 2 ' 
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The formulae (34) and (40) are changed to 



(8-5) y 6 = -4[t + l-^)u p -4y 4 , y 6 (0) = 1, 



(86) yi2 = -4^ + l-^u,-4y 10 -4, y 12 (0) = 0, 

respectively. 

Summing up the results we state, that similarly to Theorems 6 and 7, the problem of finding 
the best possible M for which the Pick function Pm locally maximizes Re a 4 in S(M) is reduced 
to the solution of the equations 



(87) y 3 (l - 1/M) = 0, 
or 

(88) 2/ 3 (l - 1/M)y 9 (l - 1/M) - y 2 15 (l - 1/M) = 0, 

where yi(t), . . . , yis(t) are solutions to the Cauchy problem for the differential equations (29-33), 
(85), (35-39), (86), and (41-43) with u p and u q given by (83) and (84). 

Theorem 12. The Pick function Pm locally maximizes Re a 4 in S(M) if M > M and does 
not give a local maximum to Re a 4 if M < Mq, where Mq is the maximum of the roots of the 
equations (87) and (88), where y x {t), . . . , y 15 (t) are solutions of (29-33), (85), (35-39), (86), 
(41-43), and (83-84). 

Numerical methods applied to (87) and (88) show that Mq ~ 22.9569... 
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